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We develop a detailed microscopic analysis of electron 
transport in normal diffusive conductors in the presence of 
proximity induced superconducting correlation. We calcu- 
lated the linear conductance of the system, the profile of the 
electric field and the densities of states. In the case of trans- 
parent metallic boundaries the temperature dependent con- 
ductance has a non-monotoneous "reentrant" structure. We 
argue that this behavior is due to nonequilibrium effects oc- 
curing in the normal metal in the presence of both supercon- 
ducting correlations and the electric field there. Low trans- 
parent tunnel barriers suppress the nonequilibrium effects and 
destroy the reentrant behavior of the conductance. If the wire 
contains a loop, the conductance shows Aharonov-Bohm os- 
cillations with the period <E>o = h/2e as a function of the 
magnetic flux $ inside the loop. The amplitude of these os- 
cillations also demonstrates the reentrant behavior vanishing 
at T — and decaying as 1/T at relatively large temperatures. 
The latter behavior is due to low energy correlated electrons 
which penetrate deep into the normal metal and "feel" the 
effect of the magnetic flux $. We point out that the density 
of states and thus the "strengh" of the proximity effect can 
be tuned by the value of the flux inside the loop. Our results 
are fully consistent with recent experimental findings. 



I. INTRODUCTION 

Recent progress in nanolithographic technology re- 
vived the interest to both experimental and theoreti- 
cal investigation of electron transport in various meso- 
scopic proximity systems consisting of superconducting 
and normal metallic layers. In such systems the Cooper 
pair wave function of a superconductor penetrates into 
a normal metal at a_jdistance which increases with de- 
creasing temperatureo. At sufficiently low temperatures 
this distance becomes large and the whole normal metal 
may acquire superconducting properties. Although this 
phenomenon has been already understood more than 
thirty years ago and intensively investigated during past 
decades, recently novel physical features of metallic prox- 
imity systemsJmye been discoveredcTcl and studied the- 
oretically (seeElEiJ and further references therein). 

In this paper we study the influence of the proxim- 
ity effect on transport properties of a diffusive conductor 
in the limit of relatively low temperatures and voltages. 
We will assume that this conductor is brought in a direct 
contact to a superconducting reservoir which serves as an 
effective injector of Cooper pairs into a normal metal. We 
will show that if the system contains no tunnel barriers 



there are two different physical regimes which determine 
the system conductance in different temperature inter- 
vals. It is well known that proximity induced supercon- 
ducting correlation between electrons in a diffusive nor- 
mal metal survives at a distance of order £jy ~ \JVjT , 
where T> = VFhmp/3 is the diffusion coefficient. As T 
is lowered the proximity induced superconductivity ex- 
pands into the normal metal and, consequently, the "nor- 
mally conducting" part of the system effectively shrinks 
in size. This effect results in increasing of the conduc- 
tance of a normal metal. At sufficiently low temperature 
the length £/v becomes of order of the size of the nor- 
mal layer and the system behavior becomes sensitive to 
a physical choice of the boundary condition at the edge 
of the normal wire opposite to that attached to a super- 
conductor. 

One possible choice of this boundary condition corre- 
sponds to the assumption that a nontransparent barrier 
is present at the edge of this wire. Then electrons can- 
not diffuse out of the wire, the proximity induced super- 
conducting correlation survives everywhere in the system 
and a real gau-in the quasiparticle spectrum develops 
in the N-metalEl The value of this gap is of order of 
e g ~ min(A, V/L 2 ), A is the bulk superconducting gap 
and L is the length of the normal wire. 

Another possible situation corresponds to the presence 
of a big normal reservoir N' directly attached to the N- 
wire by means of a highly transparent contact. In this 
case even at very low T the proximity induced Cooper 
pair amplitude is essentially nonhomogeneous in the N- 
metal. Indeed, close to a superconductor this amplitude 
is large, whereas in the vicinity of a normal reservoir it 
is essentially suppressed. Thus, strictly speaking, the 
whole N-wire cannot be characterized by the real gap in 
its quasiparticle spectrum. In the absense of a poten- 
tial barrier between N and N' this gap is obviously equal 
to zero at the NN' interface and - as will be demon- 
strated - everywhere in the normal metal. Nevertheless, 
it turns out that the density of states in the N-metal 
shows a soft pseudogap which is again of the order of 
Eg. In other words, the spatially averaged normalized 
density of states Nn(e) in the N-wire at small e < Ed 
is smaller than its normal state value iVjv < N(0) but 
always remains nonzero. It increases with increasing s 
and reaches the value Nn = N(0) at e > Ed- This is 
the key point for understanding the low temperature be- 
havior of the conductance of our system. As the tem- 
perature increases from zero, higher and higher values of 
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e contribute to the current and the system conductance 

- due to the increase of iVjv with e - increases with T. 
This regime takes place until the temperature reaches the 
value T ~ Ed where the crossover to a high temperature 
behavior takes place. Note that similar behavior of the 
normal metal conductance in the presence of proximity- 
induced superconductivity has been recently found by 
Nazarov and StoofEl 

An interesting feature of the system without tunnel 
barriers is that at T = its conductance exactly co- 
incides with that of a normal metal with no proximity 
effects. This result has been first obtained by Arte- 
menko, Volkov and Zaitsevo for the case of a normal- 
superconducting constriction. Although it is already 
around for many years the physical meaning of this result 

- if any - still needs to be understood. At the first sight 
the linear conductance of the system at T = should 
be smaller than in the normal state because of the pres- 
ence of the (pseudo)gap in the normal DOS Nn at low 
energies. Why is this not the case? 

In order to answer this question we should recall the 
well known fact that in the presence of noncquilibrium ef- 
fects the current flowing in a superconductor depends not 
only on the normal DOS but is characterized by a set of 
generalized DOS23. Our problem is just a particular ex- 
ample of a nonequilibrium superconductor: on one hand 
superconducting correlation penetrates into the normal 
metal and the Cooper pair amplitude is nonzero there, on 
the other hand in the absence of low transparent tunnel 
barriers the electric field also penetrates into the N-metal 
and drives the quasiparticle distribution function out of 
equilibrium. We will argue that in this situation one of 
the generalized DOS (below we define it as Ns(e)) - that 
is nonzero in the N-layer due to the presence of proxim- 
ity induced superconducting correlation at low energies - 
plays an important role and also contributes to the sys- 
tem conductance. In other words, in the presence of the 
electric field inside the system both uncorrelated and cor- 
related electrons contribute to a dissipative current. This 
is the reason why in the presence of proximity induced su- 
perconductivity the system conductance is never smaller 
than its normal state value although the normal DOS 
iVjv(e) < N(0) at low energiest3. 

We would like to emphasize that the situation is en- 
tirely different in the presence of low transparent tunnel 
barriers. Provided their resistances are much larger than 
that of the N-metal the whole voltage drop takes place 
at these barriers and the electric field does not penetrate 
into the N-layer. In this case only uncorrelated electrons 
contribute to the dissipative current and therefore only 
the normal DOS Nn matters. As a result the tempera- 
ture dependence of the system conductance changes. We 
will demonstrate that with lowering the barrier trans- 
parency the crossover takes place to the effective conduc- 
tance decreasing monotonously with T, characteristic for 
two serial NIS' tunnel junctions (S' is now the diffusive 
normal conductor with the proximity- induced gap). 

Note that both types of the behavior, namely reen- 



trant and monotonously decreasing withJT conductance 
have been observed in the experimentsuHrti Furthermore, 
we would like to point out that both densities of states 
iVjv(e) and Ns{e) can be investigated in one experiment. 
We will come back to this point further below. 

When the system contains a mesoscopic loop of a nor- 
mal metal, the conductance as a function of the magnetic 
flux through the loop shows oscillations with the period 
$0 = h/2e (superconducting flux quantum). Although 
the Cooper pair amplitude (and thus the supercurrent) 
in the ring is exponentially small at T > e,;, the am- 
plitude of these oscillations decays only as cx 1/T. This 
again illustrates an important difference in the behavior 
of kinetic (conductance) and thermodynamic (supercur- 
rent) quantities. Below we will argue that in the systems 
considered here even at large T 3> £d the behavior of the 
first quantity is dominated by correlated low energy elec- 
trons with e < ed penetrating far into the normal metal 
whereas the contribution of electrons with e ~ T to the 
second one is only important. Again the presence of the 
electric field inside the N-metal is crucially important 
for this effect. At low T the oscillation amplitude again 
shows the reentrant behavior and vanishes in the limit 
T — > as T 2 . A clear experimental evidence for a 1/T 
decay of the conductance oscillations has been recently 
reported in Ref. o. 

Finally we point out that making use of the the geom- 
etry with a metallic loop one can easily tune the den- 
sities of states of the system by applying a magnetic 
flux inside this loop. We will show that e.g. for 
3? = ®o/2 the proximity effect in the normal region "af- 
ter" the loop is completely suppressed and the normal 
DOS N N (e) = N(Q) there. This effect can be investigated 
experimentally and used for further studies of proximity 
induced superconductivity in normal metullic structures. 

The structure of our paper is as follows. In Section II 
we briefly describe the general kinetic approach based on 
quasiclassical Green functions in the Keldysh technique 
and define the physical quantities of interest. Then a de- 
tailed analysis of this quantities (conductance - Section 
III, DOS - Section IV, elecric field - Section V) will be 
presented. Sections VI and VII are devoted to the ex- 
tension of our analysis to the proximity systems contain- 
ing mesoscopic normal metal loops with a magnetic flux. 
The main results of the present paper are summarized in 
Section VIII. Further details related to different geomet- 
ric realizations of the proximity systems with loops are 
presented in Appendix. 

II. KINETIC ANALYSIS 

A. General formalism 

Let us consider a quasi-one-dimensional normal con- 
ductor of a length 2L with a superconducting strip of a 
thickness 2d s attached to a normal metal on the top of 
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it and two normal reservoirs attached to its edges (see 
fig.l). The length L is assumed to be much larger than 
the elastic mean free path li mp but much shorter than the 
inelastic one. This geometrical realization has a. direct re- 
lation to that investigated in the experimentsu'Btl. Two 
big normal reservoirs N' are assumed to be in thermody- 
namic equilibrium at the potentials V and respectively. 
In contrast to the case of a ballistic constrictionEZio the 
potential drop within the system is distributed between 
the interfaces and the conductor itself. The general ap- 
proach to calculate-tic conductance of these structures 
was developed inullj'113. In what follows we shall apply 
this method to analyse the temperature dependence of 
the NS proximity structure of fig. |[ 
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FIG. 1. The experimental system under consideration 

Such an experimental realization allows to prepare a 
structure without effective tunnel barriers in the direc- 
tion of the current flow. Even with "perfect" samples in 
a usual sandwich geometry, a natural barrier shows up 
due to the inevitable mismatch of Fermi velocities be- 
tween different materials. This could well be one of the 
reasons why in the previous experiments with sandwich- 
like structures, the.reentrant behavior of the conductance 
was not detectedEJ. 

The electron transport through the metallic system can 
be described by the equations for a matrix of quasiclas- 

sical Green functions G in the contac 



K 



V 

G = 



G G 

A ^ 

G 



(1) 



where G A , G R and G K are respectively the impurity- 
averaged advanced, retarded and Keldysh Green func- 
tions. These functions are in turn matrices in the Nambu 
space: 



A a A 

G =o z g + i(j y j , G = 



(G K )*and G =G f - fG 1 



Here the distribution function /= /; + a z f t , where 
fi = tanh(e/2T) and /(describes deviation from equi- 
librium. Taking advantage of the normalization con- 
dition for the normal and the anomalous Green func- 
tions (g R ) 2 — (f R ) 2 — 1 it is convenient to parametrize 
g R = coshff, f R = sinhf?, where 9 = #i +z#2 is a complex 



function. Deep in the bulk superconductor it is equal 
to 6 S = l/21n[(A + e)/(A-e)] - in/2 for e < A and 
6 S = (1/2) In [(e + A)/(e - A)] /2 for e > A (here and 
below we omit the indices R(A)). 

The current / and the electrostatic potential <f> are ex- 

v 

pressed through G as 
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vVS 



de Sp 



a\G a. G 
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/ deTig K (x,e) = / ds f t (x, e)v e (x), 
Jo Jo 



(2) 



(3) 



where v is the the density of states, v e {x) = $i(g R (x)) 
and S is the crossection area of N conductor. 

Being e™*essed in terms of the function 6(e, x) the 
equationsE3'tll for the Green functions and the distribu- 
tion function for the N-metal take a particularly simple 
form 



d x [2?(cosh 2 d x )d x f t ] =0, 



(4) 



(5) 



x is the coordinate along the N-conductor. Here we 
neglected the processes of inelastic relaxation and put 
the pair potential in the normal metal equal to zero 
Ajy = assuming the absence of electron-electron in- 
teraction there. 

Before we come to a detailed solution of the problem let 
us point out that the conclusion about the anomalous be- 
havior of the system conductance can be reached already 
from the form of eq. (||) . Indeed it is quite clear from (||) 
that the effective diffusion coefficient £>ff = 2? cosh 2 6\ 
increases in the N-regions with proximity-induced super- 
conductivity and, therefore, the electric field is partially 
expelled from these regions. This energy dependent field 
modulation is controlled by the solution for 9(e, x) and 
is directly related to the physical origin of the anoma- 
lous temperature dependence of the system conductance 
discussed below. 

The equations (f|) and (||) should be supplemented by 
the boundary conditions at the interfaces of the normal 
metal N. Assuming that the anomalous Green funjetion 
of big normal reservoirs N' is equal to zero frorrBE3 we 
obtain 

C N JBd x 9 = ±sinh 9, 
&7BCOBh0i0*/t = ±cosh^ 2 (/ t - f t (x = Q,2L)), W 

where 73 = Rb/pN^N is the interface resistance parame- 
ter, Rb is the specific resistance of the interface between 
the N-conductor and the N '-reservoirs, pn is the resistiv- 
ity of the N-metal, and £^ = ■J 'Dn /2ttT c is the temper- 
ature independent characteristic length scale in N (note 
that the coherence length in N £,n(T) = \J /2-wT is 
T-dependcnt). 
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In general we should also fix the boundary condition 
at the interface between the N-metal and the supercon- 
ductor. For the case of a perfect transparency of this 
interface and for typical thickness of the normal layer 
wn ~ \/S, Cooper pairs easily penetrate into it due to 
the proximity effect and the Green functions of the N- 
metal at relatively low energies for d < x < d + 2d s are 
equal to those of a bulk superconductor 9 = 9 S (the in- 
fluence of finite transparency of the NS-contact will be 
discussed below). In this sence the region of a normal 
metal situated directly under the superconductor can be 
also treated as a piece of a superconductor S' and the 
solution of (|j), (||) needs to be found only for < x < d 
(without loss of generality we will stick to a symmetric 
configuration). 

Cooper pairs penetrate into the normal conductor also 
in the case of a not perfectly transparent NS interface. 
As it is demonstrated below, the energy gap is induced 
in S' region in this case. As a result, for a sufficiently 
long N-wire, which is only considered here, the presence 
of the barrier at the NS interface will not influence the 
results derived for the system conductance. 



B. Physical quantities of interest 



Bipceeding along the same lines as it has been done in 
ref.113 we arrive at the final expression for the current 



I = — / de 
2RJo 



tanh 



e + eV 
2T 



tanh 



e-eV 
2T 



D{e), 
(7) 



where rQ{e) defines the effective transparency of the 
systerr£3 



D(e) = 



1 + 7" 



1 l A f L 

cosh 0i(x=O,e) cos 9 2 (x=0,e) L JO 



J Q dxsech 9i(x, e) 



(8) 



R = Rb + Rn and r = Rb/RN = 1b£,*n/L , Rpj is the 
resistance of the N-metal. 

Let us consider the case of a sufficiently long normal 
conductor d 2 ^> T) / A. Then at low temperatures T <C A 
the interesting energy interval is restricted to e -C A. 
For such values of e the contribution of the S"-part of the 
normal conductor shows no structure and can be easily 
taken into account with the aid of obvious relations 

I dxsech 2 ^i(a;,e) = / dx sech 2 #i(a;, e) + d s sech 2 Sj i 
Jo Jo 

(9) 

and sech 2 6* s .i = ( 1 — ^\ (no barrier at the NS interface) 
or sech 2 # s i = (1 — ) (the barrier is present at the NS 
interface). Due to this reason we will discuss only the 



properties of the iV-part (0 < x < d). For the sake of 
completeness we will also demonstrate the effect of finite 
d s in the end of our calculation. 

For the differential conductance of the iV-part < x < 
d normalized to its normal ("non-proximity") value in 
the zero bias limit eq. (0) yields 



Gn — 



Rdl 
~dV 



v=o 



1 

2T 



deD(e)sech 2 (e/2T). (10) 



Analogously the normalized zero-bias electrostatic poten- 
tial distribution reads 

Mx) = Hm ^ = J^J™ deD(e)v E (x)se^(e/2T) x 



x / dx'sech 2 (6»i(a; / )) 



(11) 



The normal density of states is given by the nor- 
mal Green's function via the standard relation v t (x) = 
7V(0)9fi(<? e (a;)) which enters into the conductance in the 
form cosh 9\ — (5Rg) 2 + (3/) 2 together with a 'corre- 
lation DOS' 77,(2;) = N(0)$s(f £ (x)). The importance of 
the latter quantity for understanding the effects discussed 
here has been already pointed out in the Introduction. 
We will discuss the features of these local densities as 
well as the averaged ones: 



Nn{c) — J dxv e (x) 
N s {e) = / dxrjAx) 



As it has been already mentioned the "correlation 
DOS" 7/ belongs to the set of generalized densities of 
states familiar from the-standard theory of nonequilib- 
rium superconductivityCJO. It reflects the presence of 
superconducting correlations at low energies. E.g. in a 
BCS superconductor this function reads r\ = ^ a _~^ , 
In our case this function is not only energy- but also 
space-dependent due to the fact that the proximity in- 
duced superconducting correlation decays inside the nor- 
mal metal. But the physical meaning of it remains the 
same as. in standard noncquilibrium superconductivity 
theoryc3: 77 plays a role whenever the quasiparticle dis- 
tribution function of a superconductor is driven out of 



known problems 
It happens also 



equilibrium. It happens e.g. in the wi 
of charge relaxationEil and imbalance 
here due to a simultaneous presence of the electric field 
and the proximity induced superconducting correlation 
in the normal metal. 



C. Influence of finite barrier transparency at the top 
NS-interface 

Let us consider the effect of a tunnel barrier at the NS 
interface in more details. 
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Under the assumption that the N-wire thickness is 
small Wn *C £jv the equation for 9 in the region < 
x < 2d a underneath the superconducting terminal can 
be derived by the method of Ref.EJ: 

Pd 2 6> + 2iesinh6> + Acos(9 = (12) 

where the effective order order parameter A = 
sin6> s /7^ s , the effective energy e = e + cos 8 s /~fg S , 
7s 5 — {fib I ' Pn){wn /£n) * s the interface transparency 
parameter. Here 9 S is the solution in S which is set 
equal to the bulk value tan _1 (iA/e), a good approxi- 
mation for thin N film wn <C £ n . With these substitu- 
tions the equation ([l2]) in the N film has the form similar 
to that in a superconductor. This equation is valid for 
7s S > { w n/£,n) 2 i- e - f° r sufficiently small transparency 
of the NS interface: (D) < In/wn- 

As follows from Eq. (|l2|) , superconducting properties 
of the N-layer are described in terms of the energy- 
dependent coherence length 

1 /2 

6v(e) = {hD N /2 [ lB 2 -e 2 - 2i lB \ cos 9 S ] 1/2 } , 

(13) 

which determines an exponential decay of Nn(x, e) with 
x . 

The expression for £/v(e) has a pole at the gap energy 
e = A 9 jVi which signals the decay of quasiparticles en- 
tering N at e < A gN . At high energy e > ttT c /^ s 
the well known result £z\r(e) = {^Dn /2e) x / 2 is repro- 
duced, whereas at low energies e <C irT c /jg S one obtains 
^ = {hD N ^ s /2nT c ) 1 / 2 . Thus the effective length 
scale in N increases with the decrease of the SN inter- 
face transparency. 

It is straightforward to calculate the gap energy 
A g N assuming the "rigid" boundary conditions 9 S — 
tan -1 (iA/e) either from the pole of £/v(e) or, making use 
of the solution 9 = tanh - [smh9 s /(cosh9 s — ie^^ s )] 
and calculating the quasiparticle density of states 
N N (e) = Recosh9. Subsituting the expression for 9 S 
into this solution one arrives the equation for the energy 
gap 

t 3 + 2Ct 2 + (C 2 - l)t - 2C = 0, (14) 

where t = yj\ - (A gA r/A) 2 and C = t:T c /j% s A. The 
general solution is rather cumbersome, therefore here we 
present only its asymptotic forms. 

The gap is given by A gN /A = 1 - 2( 7 f s A/ttT c ) 2 for 
large transparency of the NS interface, j B s A/ttT c <C 1, 
and by A„ N = A/(l + j^ s A/nT c ) for small trans- 
parency, jg S A/TrT c > 1. In the latter case (the McMil- 
lan limit), the expression for the gap may be written 
as AgN — irT c /^g S = Hvfn {D} /4u>jv. The gap origi- 
nates from the finite average lifetime tn = 2wn/ (D) vfn 
for quasiparticles in the N layer with respect to Andreev 



scattering from S, since a contribution of gapless quasi- 
particle trajectories parallel to the NS interface is elim- 
inated in the diffusive regime. The dependence of A g N 
on 7g in the whole range of transparencies is presented 
in Figj2| Note that under a substitution wn i^jr/kp the 
above expressions reproduce the result of RefH for a gap 
induced in a 2D electron gas in contact to a superconduc- 
tor. We also note that in the case of two superconducting 
terminals attached to N the gap aquires a phase factor 
A g jv cos(f/2)A gN , ip being the phase difference be- 
tween the terminals. 

At the subgap energies e < A g N the solution in N is 
9 ~ l/21n[(A g jv + £)/(A 9 jv — e)] — in/ 2 whereas above 
the gap, at e > A gN , 9 ~ 1/2 In [(e + A gN )/{e - A gN )}, 
in complete analogy with the solutions in a bulk super- 
conductor but with reduced gap. Therefore the results 
for conductance of N depend on the relation between the 
gap A g and the Thouless energy Ed = D/d 2 . Since for 
the case of a sufficiently long N-wire, which is only con- 
sidered in this paper, the condition A g N ^ £d is satisfied, 
the presence of the barrier at the NS interface does not 
influence our results for the system conductance. 




FIG. 2. Reduced gap on the normal side of the 
iVS-boundary in the presence of barriers 

III. CONDUCTANCE 

A. Perfectly transparent boundaries 

The analysis of the problem can be significantly sim- 
plified in the case of perfectly transparent interfaces 
(7b = 0). In this case the boundary conditions are 

6»(0) = (15) 
0{d) = 9 S (16) 

for the contact to the normal and the superconducting 
reservoir respectively. The effective transparency of the 
N-part then reads 

D(e)=(~[ dxsech 2 (9 1 (x))\ . (17) 
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As it was already pointed out for relatively long normal 
conductors and at low T only the energies e <C A give 
an important contribution to the conductance. In this 
case the typical energy scale is defined by the Thouless 
energy = T>/d 2 -C A, A g iy. For these energies we can 
set 9s = —iir/2. Let us first put T = 0. Then the thermal 
distribution factor sech 2 (e/2T)/(2T) reduces to a delta 
function and we have 



St(fc)(l-z) 



G N (T = 0) = D(Q), 



(18) 



i.e. we only need the solution of (Q) with boundary con- 
ditions ( |l5| ) at e — 0, which is 9 = —i^x. This does not 
depend on V, so the correlations are destroyed by the 
influence of the boundary conditions but not by thermal 
excitation or by impurity scattering. From here, we can 
calculate the conductance 



G N {T = 0) = 1, 



(19) 



i.e. at T = the system conductance exactly coincides 
with its normal state value (cf.cJ'Ea. This result, however, 
by no means implies the destruction of the proximity 
induced superconductivity in the N-layer. Later on, we 
will demonstrate, that the DOS and the electrical field 
are completely different from their values in the normal 
state and in fact only due the additional contribution of 
correlated electrons the system conductance does not fall 
below its normal state value. 

In the case T -C £d we can calculate 9 perturbatively. 
From Vd 2 x 9 = -2ie sinh6» (x) and @ we get 

8 e , , ._ ... 7r _ 

9 = 7T — [x — sm(a;7r/2) — %—x. 

7T s d 2 

Keeping only leading order terms in ~, we get 



0JV 



1 d 



(20) 



where A — (I — ^) « 0.049 is a universal con- 
stant. This means, that for low temperatures Gn(T) 
grows quadratically on the scale of Ed and approaches 
the crossover towards the high temperature regime dis- 
cussed below. 

In the limit T » e d (where we still have T -C A), 
the contribution of the low energy components to the 
thermally weighted integral for G_n(T) is oc 1/T as we 
will see below and can therefore be neglected. We only 
have to take into account the solutions of (||) for energies 
e ^s> Cd- It is well known (see e.g.113), that for this energy 
range the solution of (^) together with (jl5|) reads 



tanh(0(x)/4) = tanh — e 



fe(x-l) 



(21) 



where k = dyl — 2ie/T>. By using obvious substitutions 
and multiple-argument relations for hyperbolic functions, 
we arrive at the following identity: 



dx sech 




q{y) dy 



(i 



q(y)) 2 

(22) 



where q(y) = 4(3 + 2V2>- 2 ^ sin 2 yj [eT^ + 3 + 2^2) . 

For calculating D{e) we can, as the integrand becomes 
exponentially small for y > §t(k) S> 1, take the upper 
bound to infinity, such that it becomes a universal con- 
stant. From there we can calculate the conductance in 
this limit 



G N (T) 




(23) 



where again B = 0.42 is a universal constant. 

These results has a simple physical interpretation. Su- 
perconductivity penetrates into the normal part up to 

£w = J 2^!p, whereas the rest stays normal, so the total 
voltage drops over a reduced distance d — £n- Thus the 
resistance of the structure is reduced according to the 
Ohm law. In terms of the conductance, this means 



G N = l + B 



(24) 



which is equivalent to 




FIG. 3. Local effective diffusion constant 

Let us point out, that at both edges of the N- metal the 
local effective diffusion constant £> c ff = cosh 2 9\D is not 
enhanced (see Fig. ||) in comparison to its normal state 
value, because either the Cooper pair amplitude (at the 
NN' boundary) or the electric field (at the NS boundary) 
is equal to zero due to the imposed boundary conditions. 
Inside the N-metal the value 2? e g- becomes higher due to 
nonequilibrium effects in the presence of superconducting 
correlations (77 ^ 0). This effect is small at very low 
energies and becomes more pronounced at e ~ Ed- 

For temperatures comparable to Ed the problem was 
treated numerically. The results show an excellent agree- 
ment with our analytical expressions obtained in the cor- 
responding limits. 
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FIG. 4. Conductance in the case of transparent barriers 

The numerical results (see fig. EJ) confirm, that for 
Ed <C A the universal scaling with ^ is excellently ful- 
filled, the conductance peak with the height of about 9% 
takes place at T w 5£^ (cfJ23). This peak becomes smaller 
if we take into account the influence of finite ds keeping 
d fixed (fig. |^) The qualitative features, however, remain 
the same. 
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FIG. 5. Conductance normalized to the total length of the 
normal wire, rj := ds/d, td = 10 A 



B. Tunnel barriers 

Let us now assume that a tunnel barrier is present 
at the N'-N interface. If one lowers the transparency 
of this barrier the crossover takes place to the behav- 
ior demonstrating monotonously decreasing conductance 
with T (Fig. ||), which is typical for two serial NIS tun- 
nel junctions. Fig. |^ demonstrates the crossover with 
increasing r — 7s£^/d. Inset shows the Arrenius plot 
for the case of 7s^/(i ^> 1 which illustrates the acti- 
vated tunncl-likc behavior. 
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FIG. 6. Conductance in the presence of tunneling barriers 

Formally this is due to the term r/ cosh #1(2; = 
0, e) cos Qi(x = 0, e) in the denominator of Eq.(||), i.e. the 
contribution of the barrier at the N'-N interface. In the 
small transparency limit r 1 the expression (||) reduces 
to the standard tunnel formula. The physical reason for 
this behavior is transparent. For r <C 1 the presence of 
a tunnel barrier is not important, the electric field pene- 
trates inside the normal metal and we come back to the 
picture discussed above for perfectly transparent bound- 
aries in which case both normal and correlation DOS play 
a significant role. If, however, the resistance of a tunnel 
barrier dominates over the Drude resistance of the nor- 
mal metal r>l, the whole voltage drop is concentrated 
at the barrier, nonequilibrium effects in the N-metal are 
absent and therefore only the normal density of states 
enters into the system conductance. 

An additional effect is that a real gap instead of a soft 
pseudogap develops in the case low transparent tunnel 
barriers. The crossover between these two regimes is dis- 
cussed in more detail below. 

Note that both types of behavior, namely nonmonoto- 
neous and monotonously decreasing with T conductance 
have been observed in the experiments!! 
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FIG. 7. Averaged normal and correlation DOS 




IV. DENSITY OF STATES 



A. Averaged density 

From our approximative solutions of the preceeding 
sections, the densities of states can be easily calcu- 
lated. For e = we have N N = N s = N(Q)%. 
At low energies e « q there are quadratic correc- 

tions: N N/S = N(0) U ± A 1/2 with A x = 



16 



'1 + ^1 « 0.198. 



« 0.0396 and A 2 
For high energies, the densities approach their nor- 
mal values, again with square-root corrections Njy = 
N(0) (1 - Biypf) and N s = B 2y ^f with B x « 0.321 
and B 2 » 0.75. ' 

Together with our numerical data (see Fig. 0), this 
demonstrates the presence of a soft pseudogap in the den- 
sity of states below the energy Similar results have 
also been discussed inE3. 



FIG. 8. Local DOS for different N-N'-boundaries. 
Transparent, Bottom: Non-transparent 



Top: 



B. Spatial dependence and tunneling experiments 

It is also interesting to investigate the spatial depen- 
dence of the DOS in the normal layer. Fig. || show lo- 
cal normal DOS Nn calculated for perfectly transparent 
(7s = 0) and nontransparent (jb = °o) NN' interfaces, 
respectively, at different distances from the NS boundary: 
x/d= 0.1, 0.2, 1. The difference between these two 
cases is quite obvious: whereas for — the normal 
DOS at low energies is always finite, becoming larger at 
larger values of x, for 75 = cxd a real gap in the density 
of states clearly shows up at all energieSj-Similar results 
have been recently discussed in RefsH3ili. The overall 
behavior of the local correlation DOS at each value of x 
is similar to its average value. 

It is important to emphasize that both v{x) and rj(x) 
are measurable quantities and can be directly probed 
in experiments. Recently the spatial and energy depen- 
dence of the normal DOS has been studied in tunneling 
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experiments^. The data@ show a qualitative agreement 
with theoretical predictions. The results obtained here 
suggest that much better agreement can be achieved if 
one takes into account smearing of the proximity induced 
gap in the normal metal due to the diffusion of normal 
electrons from the external circuit (which plays the role 
of the N' reservoir) through the NN' boundary. For non- 
transparent NN' boundaries (75 = 00) this process can 
be neglected and a real gap develops in-ljhe N-metal (Fig. 
|§| (b)). As no such gap was found irO we believe that 
diffusion of normal excitations into the N-layer from the 
external circuit should play an important role in these 
experiments. In other words, the experimental situation 
appears to be closer to that described by the boundary 
condition 7s = with a soft pseudogap (Fig. || (a)) than 
to the case 75 = 00 (see e.g. Ref. |4l| ). The dependence 
of this effect on the size of the N-layer is depicted in Fig. 



dl_ 

dV 



2tt 2 T 2 



(26) 



v=o 



These our results demonstrate that the depairing effect of 
the N'-reservoir needs to be taking into account on cq.ua.1 
footing with pairbreaking due to inelastic scatteringOEj. 

Let us also point out that one can also extract infor- 
mation about the correlation DOS by making two kinds 
of measurements with the same sample. Indeed, by mea- 
suring the conductance of the system (or a part of it) 
with no tunnel barriers one obtains information about 
the combination of Nn and Ns entering the expression 
for the system conductance G, whereas performing the 
tunnel experiments^] one probes only the normal DOS 
Nn. Then the correlation DOS can be easily recovered. 



V. ELECTRIC FIELD AND CHARGE 




FIG. 9. Size effect on the local normal DOS. Here, the 
density of states at a fixed distance x = 0.5 from the 
iVS-boundary is plotted for different values of the total length 
L of the N-part. 

Making use of the Usadel equation one can easily re- 
cover simple analytic expressions for the density of states 
at a distance £0 away from the NS-boundary. For a N- 
wire of the total length d at e <C we obtain 



N N (e,x ) 



e 2 

where a and (3 describe the size effect 



(25) 



d 



■ f 71 

a = sm — 
V2 

/3 = ^ a ( 1 "T 



'TTXq 

S, 2T 



In this section we shall discuss only the case of perfectly 
transparent interfaces. 

^From our solutions we can calculate the electric field 
and the charge by using ( |Tl| ) and the Poisson's equation. 
The field shows essentially non-monotoneous behavior. 
At T = we have E(x) = cos(z7r/2)-§ (x-1) sin(x7r/2). 
At high temperatures T ^> ed, the field is constant E = 1 
far from the superconductor where no correlation remains 
(I — x ^ £ e ) and it changes linearly near the supercon- 
ductor: E{x) — -64(1 — %)\J~^2 with _B 4 w 2.59, however, 

it still overshoots in between these regimes (see Fig. [To| ). 

We see that close to the superconductor the electric 
field monotoneously decreases with temperature as su- 
perconductivity becomes stronger there. Further from 
the NS boundary the field shows a complicated behav- 
ior overshooting the normal state value (the total volt- 
age drop is fixed!) in the region where superconducting 
correlation starts decaying either due to thermal effects 
(high T) or due to the presence of a normal reservoir 
(low T) . The local resistivity is maximally lowered there 
and the layer of polarization charges is formed (see Fig. 
|To| ). These results emphasize again the importance of 
nonequilibrium effects for understanding the behavior of 
the system conductance. 



Thus for xq <C d the normal DOS at zero energy and 
x — xq is proportional to X/d. Neglecting the charging 
effects (which in principle can also be importantcS) for 
the differential conductance of the tunneling probe we 
find 
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If the wire was a real superconductor, the magnetic flux 
would induce a supercurrent into the ring. As a function 
of $, this current has a period of the superconducting flux 
quantum $o = h/2e. 

To describe these type of systems, our kinetic scheme 
has to be extended in several points. 

We define the Green's functions in the loop as 

G = coshw e F = smhu e e liPe <p e = Xe - 2e / dl A(l) 

Jo 

where the integration goes along the loop. In the pres- 
ence of a vector potential, we have to introduce gauge 
independent derivatives 

V — > V - 2ieA. 
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(a) Distribution of the electrical field 
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(b) Distribution of the electrical charge 
FIG. 10. Electrostatics within the wire 



VI. EXTENSION TO SYSTEMS CONTAINING A 
LOOP 



Recently, the properties of proximity wires 
ing a loop have attracted much experimenta 
theoreticalcjO interest. 



t ain- 
and 




FIG. 11. The system under consideration. S and N are 
respectively superconducting and normal reservoirs. The wire 
is made of normal conducting material. 



This means, that instead of solving the Usadel equation 
with a vector potential, we can perform a gauge trans- 
formation and map onto a system without magnetic field 
having phase x instead of <j). As the definition of the 
Green's functions has to be unique everywhere in the 
loop, we have lim K _ > o+ vix) — lim^^o- ( mod 27r) 



lim X (x) 



lim x( x ) 



2e$ 



( mod 2ir) 



after gauge transformation. Here, $ is the magnetic flux 
in the ring. 

This mapping shows, that the magnetic field induces a 
supercurrent j~ (screening current) into the system. We 
want to neglect any conversion between this supercur- 
rent and the dissipative current, so both are conserved 
seperately. This allows the application of the kinetic 
scheme which h&s been developed for systems without 
phase gradientt2l but can be generalized to any system 
where the dissipative current is copspsved. 

The Usadel equation then readscaO 



V— 

dx 2 c 



0- • V, ^ V ( d *e 

-2iesinnu e H — — 

2 \ dx 



sinh 2u f 



(27) 



and has to be solved together with the equation for the 
conservation of the supercurrent 



TJe =0 3? = |smhu £ | — 
dx dx 



(28) 



In order to match the Green's functions at branching 
poits we use the standard continuity condition as no tun- 
nel barriers are assumed to be there. ^From the Usadel 
equation in matrix form 

DV(g e V9e) + ie[T z ,ge\ =0 
follows for any branching point (see also@) 

d 
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where the sum runs over matching branches, -J^ denotes 
the derivative in the direction of branch i and A t is the 
cross-section area of branch i. Using our definitions, we 
get 



du 



9Xe 

dxi 







(29) 



These conditions are equivalent to current conservation, 
so this is a "Green's functions Kirchhoff law" . For N = 1 
Zaitsev's boundary conditioned for a Normal- Vacuum 
boundary is reproduced, N — 2 is equivalent to the triv- 
ial statement, that the Green's functions' derivatives are 
continuous within a branch. 

For the calculation of the total transparency D = 1/m, 
we can use the fact, that the to.; fulfill Ohm's law just by 
their definition: 



d\TYt\ + d^rri^ + ( 



d.2m2 d 3 m 3 



(d 1 + (l/d 2 + l/d 3 )- 1 +d i ) 



VII. MAGNETORESISTANCE OSCILLATIONS 



The equations (|27]) and (|28|) together with boundary 
conditions (^) and branching conditions (^9|) have been 
solved numerically and also analitically in some limit- 
ing cases. For numerics, the problem was mapped onto 
a simpler boundary value problem without any fitting 
point. As the systempef equations is unstable, we used 
the relaxation methodEJ instead of shooting. 

For convenience, we have chosen d\ = d 2 = da = d± 
and Ai — 2A 2 = 2 A3 = A4, which simplifies the condi- 
tions |2^. The effect of geometry on the conductance oscil- 
lations will be discussed in the appendix. The Thouless 
energy of just one branch will be labeled as e c i = ^i- 



A. T-dependent Amplitude of /i/2e-Oscillations 

For T = 0, only quasiparticles with the energy e = 
contribute to the conductance. From ^ and g we can 
conclude, that w e =o is a purely imaginary function, so the 
total conductance of the system is equal to its normal 
state value, being independent of $. In other word s, 
there exist no conductance oscillations at T — (cf.tScZl) . 

At nonzero temperatures the system conductance de- 
pends on the magnetic flux inside the loop with the pe- 
riod equal to the flux quantum <& (see fig. |l^). With 
the aid of simple analytic arguments (see Appendix) one 
can conclude that at low temperatures the amplitude of 
the conductance oscillations increases as T 2 (see fig. 13). 
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FIG. 12. /i/2e-periodic structure of the conductance, nor- 
malized to the zero-field conductance 

In order to establish the temperature dependence of 
this amplitude at higher T 3> td it is convenient to make 
use of the fact that for electrons with sufficiently large 
energies e > td superconducting correlation is destroyed 
already before they reach the loop. Thus at such energies 
the transparency of the whole structure D(e) should be 
insensitive to the particular value of the flux inside the 
loop. In other words, calculating the flux-dependent part 
of the system conductance we can take into account only 
the contribution of low energy quasiparticles which re- 
main correlated in the loop area. E.g. for the amplitude 
of the conductance oscillations we get 



AG(T) = G h/ UT) - G (T) 
1 



2T 



2T 



-AD, 



de {D h/ie (e) - D a (e))sech'(e/2T) 

de (D h/4e (e) - A)(e))sech 2 ( e /2T) 

(30) 



where e c is the cutoff parameter of order e^, and AD a 
is constant. 
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FIG. 13. Oscillation amplitude 
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(a) Transparency (b) Conductance 

FIG. 14. Temperature dependence of tranport properties 
at different fluxes 

The results of our numerical analysis fully support the 
simple analytic arguments presented above. The system 
transparency D(e) is depicted in Fig. |l4| (a) for different 
values of the flux $. The value D(e) depends on $ only 
at low energies, whereas for e > e<j all curves merge. In 
accordance to our simple estimate ([30j) this leads to a 
1/T decay of the oscillation amplitude AG at large T 
(see Figs. |l| (b) and Also the T 2 behavior of AG 
in the low temperature limit is recovered (Fig. [l2|) . 

The 1/T behavior of AG has been also found in recent 
experimentsS'Ll. We would like to point out that a slow 
power-law decay of the conductance due to a dominating 
contribution of low energy quasiparticles just emphasizes 
the physical difference between kinetic and thermody- 
namic quantities, like supercurrent which decays expo- 
nentially with increasing T. 

B. Flux-Dependent DOS 

As it was already discussed the simultaneous presence 
of correlated electrons and the electric field in the nor- 
mal metal causes nontrivial nonequilibrium effects, the 
description of which involves two densities of states v{x) 
and r)(x). In the presence of the normal metal loop with 
the magnetic flux $ in our system there appeas a possibil- 
ity to tune both normal and correlation DOS by changing 
the value of For the system depicted in Fig. [ll] these 
densities of states can be easily calculated. As one might 
expect for the region between the superconductor and 
the loop (between the points A and B) this dependence 
is quite weak and both DOS practically coincide with 
those calculated above for a wire without the loop. On 
the other hand, in the region between the loop and the 
normal reservoir N' (between the points C and D) the 
quantities v{x) and rj(x) are very sensitive to the flux <E>. 



(a) Normal (b) Correlation 

FIG. 15. Flux-dependence of the two densities of states in 
point C of the system 

The normal and correlation DOS for the point C (Fig. 
|j"l"l) are presented in Fig. [l^. We see that with increas- 
ing the value of the magnetic flux the proximity induced 
pseudogap decreases and vanishes completely as the flux 
approaches the value $o/2. For such value of $ the prox- 
imity effect in the region "after" the loop is completely 
destroyed, the pseudogap is fully suppressed and the nor- 
mal DOS coincides with N(0) at all energies. Accord- 
ingly the correlation DOS vanishes at $ = < &o/2. Thus 
in this case the resistance of the region between the points 
C and D is equal to its normal state value at all T. 

These results demonstrate that "the strength" of the 
proximity effect in our system can be regulated by the ex- 
ternal magnetic flux. This might serve as an additional 
experimental tool for investigation of proximity induced 
superconductivity in normal metalic structures. In par- 
ticular we believe that it would be interesting to repeat 
the tunneling experimentsc3 in the presence of the loop 
with the magnetic flux. Such experiments would provide 
a direct probe of the dependence of the densities of states 
on $ (Fig. ||(a)). 

VIII. SUMMARY AND OUTLOOK 

We have used a microscopic kinetic analysis to describe 
the transport properties of superconductor-normal metal 
proximity structures. In the case of transparent inter- 
metallic boundaries we demonstrated a reentrant behav- 
ior of the system conductance with temperature. This be- 
havior was attributed to nonequilibrium effects occuring 
in the normal metal in the presence of proximity induced 
superconductivity and the electric field. We argued that 
under these conditions both uncorrelated ( "normal" ) and 
correlated ("superconducting") electrons contribute to 
the system conductance which is henceforth defined by a 
combination of two densities of states - the normal DOS 
Nn and the correlation DOS Ns- The latter is known 
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to play an important role whenever the quasiparticle dis- 
tribution function in a superconductor is driven out of 
equilibriumEj. 

We studied the energy dependence of both these quan- 
tities and demonstrated that if the normal metal is 
brought in a direct contact to a superconductor on one 
side and a big normal reservoir N' on the other side the 
normal DOS Nn shows a soft pseudogap at energies be- 
low the Thouless energy e<j- This effect is due to an 
interplay between correlated and uncorrelated electrons 
penetrating into the N-layer respectively from a super- 
conductor and a normal reservoir. If a low transparency 
tunnel barrier is present at the NN' interface the diffusion 
of normal excitations into the N-metal is suppressed, the 
influence of a superconductor prevails and a real gap in 
the density of states develops. 

Low transparent tunnel barriers also prevent the elec- 
tric field from penetration into the N-layer thus sup- 
pressing noncquilibrium effects there. We demonstrated 
that with the aid of a proper combination of the systems 
with and without tunnel barriers one can directly probe 
both energy and spacial dependencies of both densities 
of states Nn and N$ in one experiment. 

We extended our analysis to proximity systems con- 
taining the normal metal loop with the magnetic flux <E>. 
We demonstrated that the conductance of such systems 
as a function of <& oscillates with the period equal to the 
flux quantum <£>o = h/2e. The amplitude of these oscilla- 
tions AG also shows the reentrant behavior being equal 
to zero at T = 0, increasing as T 2 at T < and de- 
caying as 1/T at T > erf. We argued that even at high 
temperatures T > low energy electrons with e < 
are only responsible for the conductance oscillations lead- 
ing to the power law decay of AG at large T. We pointed 
out that the densities of states Nn and N$ can be tuned 
(decreased) by applying the magnetic flux <£>. In partic- 
ular, if the flux in the loop is equal to the half of the flux 
quantum $ = $o the proximity effect in the region "af- 
ter" the loop is completely suppressed, the normal DOS 
N N = N(0) at all energies and N s = 0. This effect 
can be also directly probed in tunneling experiments and 
used for further studies of proximity induced supercon- 
ductivity in normal metallic systems. 
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C.Bruder, W.Belzig, H.Courtois, D.Esteve, B.Pannetier, 
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van Wees. This work was supported by the Deutsche 
Forschungsgemeinschaft within the Sonderforschungs- 
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RFFI No.96-02-1956. 



containing a mesoscopic normal metal loop with the mag- 
netic flux (see Fig. |ll|). In the first three sections we will 
keep Ai = 2A^ = 2A^ = A4 for simplicity allowing differ- 
ent values for the di but restricting ourself to symmetric 
loops d2 — C?3. 



1. Low temperature behavior 

For examining the low energy range, which is dominant 
for the conductance oscillations at any temperature as 
stated above, we proceed perturbatively from the case 
e = 0, <f> = first to finite flux, then to finite e. As the 
Usadel equation is quadratic in the value u does not 
distinguish between the upper and the lower branches 
of the ring (2 or 3 in Fig. [ll]), so we will not make a 
difference in the notation. 

We start from e = 0, $ = 0, where u is purely imag- 
inary and therefore yields D = 1. For a finite value of 
the flux but e = 0, we get a purely imaginary correction 
and therefore -D^o — 1- This correction is quadratic in 
the flux as the r.h.s. of the equation is quadratic in the 
phase. The finite energy correction at zero flux is a real 
function and is quadratic in e. 

Thus proceeding perturbatively we find 5Ru oc e<i> 2 , and 
from the expansion 



u = -i-Trr + i<l> 9di,d2,d 3 ,di( x ) + thd^d^d^d^x) 
My, 

+e(t> 2 kd u d 2 ,d 3 ,dA x ) + it<t> 2 ldi,d 2 ,d 3 ,di( x ) 



with ds — d\ + d2 + d± , in the leading order we get 



<fe Jo cosh (ui(x)) I 



where the coefficients are defined as 



A= / dxbf[x) 
Jo 

fi = 4 / dxh(x)k(x). 
Jo 



APPENDIX A: A CLOSER LOOK AT THE 
H/2E-OSCILLATIONS 

Here we present further details related to the effect of 
geometry on the behavior of the proximity NS systems 



Thus at low T both the transparency and the conduc- 
tance depend quadratically on energy,and flux. Further 
analytic expressions are presented irJ13. 

It is remarkable, that for rf 3 = due to @ we have 
sinhu = at the point G. Therefore, the current conser- 
vation condition (28) can be fulfilled for j = 0, x = ±0/2 
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(different signs refer to different branches), so the phase 
gradient is zero almost everywhere and the Usadel equa- 
tion does not contain the phase any more. Thus no mag- 
netoresistance oscillations occur in this case. 



2. Cutoff energy 

Let us estimate the cutoff energy e c . 
Consider the case e 3> (-d x +d 2 — (d 1 +d 2 y 1 • ^ or ^ = 
we again have ( |2l| ) u = 4arctan f^e~ fcB:c J , so for small 

4> we can proceed perturbutavely. As the supercurrent is 
exponentially small, we can approximate the phase pro- 
file as 

d X _ 1 c k 1 (x-d 1 ~d 2 /2) 

dx 2sinhfcic?2 

so the influence of the magnetic flux is concentrated 
within a distance max{^7v,e, ^2} from the point C (see 
Fig. [n]). However, as u is exponentially small there, 
the oscillations of the transparency are exponentially sur- 
pressed, so we can estimate e c < e c i 1 +d 2 , which depends 
only on the sum of these length, but not on d\ alone (see 
also fig. ||). 



3. Cross section effect 

For the sake of simplicity above we have sticked to 
the case 2A 2 = A\ = A3. As this condition might 
not be fulfilled in real experiments it is worthwhile to 
check whether the main features of our analysis survive 
for other values of Al,2,3- in order to do that we per- 
formed calculations also for the case A\ — A2 = A3. The 
results are similar to those obtained before, showing an 
additional dip structure in the transparency at interme- 
diate energies (Fig. [l7| (a)) and a slightly deformed G in 
the same energy interval ( Fig. p"7| (b)). 
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FIG. 17. Transport quantities for the system with 
Ai — A2 = A3 — A4 and di = d? — d$ — d± displaying 
the cross-section effect. 



(a) Transparency 



(b) Amplitude of conduc- 
tance oscillations 



FIG. 16. Size effects onto 
tions. Small: di = di — 
Medium:di = = ^3 = d<t = 
d 2 = d 3 1.4. 



the conductance oscilla- 
1.2, d 2 = d 3 = 0.6, 
1, Big: di = d 4 = 0.8, 



The key feature of Fig. [Uj (b) is the strong flux- 
dependence for systems with small rings. This fact can 
be also recovered from the Usadel equation: for d 2 <C d 

As this enters quadratically 



one estimates 



d X 



dx 2<E> 

and in the end we have to integrate over the ring only 
once, the contribution of the ring is roughly oc l/da- 
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